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Tki s is the first Semi~Axuro.il status report describing research on molecular 
processes in a high temperature shock layer. 

I. Introduction 

Models 1 of the shock layer encountered by an Aeroassisted Orbital Transfer 

Vehicle require as input accurate cross sections and rate constants for the atomic 

and molecular processes that characterize the shock radiation. From the estimated 

atomic mid molecular densities^ in the shock layer and the expected residence time 

of 1 msec ^ , it can he expected that electron-ion collision processes will be 

important in the shock model. Electron capture by molecular ions followed by 

♦ — 

dissociation, e.g. 0. ♦ e -* 0 + 0 , can he expected to be of major importance 

since these processes are known to have high rates (e.g. 10 ^cm^/sec) at room 

temperature, however, there have been no experimental measurements of dissociative 
recombination (DR) at temperatures (>12000R) that are expected to characterize the 
shock layer. Indeed, even at room temperature, it is often difficult to perform 
experiments that determine the dependence of the translational energy and quantum 
yields of the product atoms on the electronic and vibrational state of the reactant 
molecular ions. This report presents ab initio quantum chemical studies of DR for 
molecular ions that are likely to be important in the atmospheric shock layer. 

II. Summary of Prior Research 

A theoretical study of DR involves the determination of potential energy 

3 ... 

curves for the molecular ion and the dissociative states of the neutral molecule. 

In order to determine dissociative recombination rates and cross sections one must 
calculate the electronic probability for electron capture given by Fermi's Golden 
Rule , 


r =2 jrf <P VlH [Q I 


where* is a mu It icon f iguration wave function and P and Q are Feshbach projection 
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Operators that project onto the tens in 9 which represent the molecular ion plus 

a free electron and the tens that represent the auto ionizing state respectively . 

4 

In the usual Feshbach projection operator formalism an eigenvalue equation is 
derived for P9 which is difficult to solve since it contains an energy dependent 
optical potential. However. I have developed a new technique in which P* can 
easily be detenined by solving the usual configuration interaction (Cl) problem. 
Vriting the total wave function as f = P9 ♦ Qt we can write the Schroedinger 
equation in matrix form as 



( 1 ) 


where Hj, p = PHP, H Qp = QHP, etc. 


Multiplying the matrices in (1) leads to: 


Hp p P * + ®i> q Q 9 = EP f and (2) 

H P 9 + H Q 9 = EQ 9. (3) 

QP QQ 

From (3) we have. 


Q*-yp*r<E - V 


(4) 


Substituting (4) into (2) leads to a matrix optical potential for P % 


(H PP * W (E - "qq” p ’ ' EPt - 


(5) 


Since E is on both sides of Eq.(5) it is difficult to solve directly for P9 . 

However P9 can be easily determined by simply diagonalizing the full H matrix in 

(1) and retaining only the coefficients of the P9 configurations. In order to 

divide by E - H in (4) and to be certain that Q does not mix into the Pi 1 that is 
vV 

determined by diagonalizing the H matrix in (1) it is necessary to first solve for 
the Q9 roots by diagonalizing H . The low energy Q9 roots are then projected out 
of the H portion of the H matrix in (1). P9 is then determined by diagonalizing 



the transformed 1 aatrix. The free electron is represented by s Rydberg orbital 
with a high principal quantum saber. Successive It's are determined for a aeries 
of increasing principal quantum numbers. The widths obtained by this procedure are 
then extrapolated to the continuum to yield the free electron capture width. 

III. Comparison of Ab Initio and Results 

During the past six months these techniques have been tested on the HO 

molecule where there are experimentally derived^ interaction matrix elements for 
2 

the B and L D repulsive states of HO. These states are expected to play an 
important role in the DR of the ground state of HO 4 . Using the approach outlined 
above for calculation of the entrance and exit channel wave functions and the 
widths, the importance of correlation of the NO 3 o and 4o orbitals has been 
investigated in the context of a full valence space Cl. These orbitals correspond 
to the atomic 2s orbitals of N and 0 at large internuclear separation. Inclusion 
of valence correlation of the 3o and 4o orbitals in the repulsive or autoionizing 
states increases the B state width and decreases the L state width by about 10Z and 
40Z respectively. In the P space entrance channel the 3 o and 4 o correlation 
decreases the B and L state widths by 12Z and 39Z respectively. This inner shell 
correlation is unimportant for the B state but very important for the L state 
width. Inner shell correlation has therefore been included in the NO width 

calculations . 

Gallusser and Dressier^ have reported width matrix elements for the lowest NO 

2 

n Rydberg states. Even though the procedure outlined above uses the high (n=7,8) 

Rydberg states for extrapolation to the continuum width it is nevertheless 

instructive to compare our calculated matrix elements for the lowest Rydberg states 

2 

to the experimentally derived results. For the B n repulsive state Hamiltonian 

2 5 

tn'Lrix element with the 3pn C D Rydberg state the calculated and experimental 

-1 -1 2 

results are 1229cm and 1382.6cm while for the 4p nK HRydberg state the results 
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are 714.2cm and 803.9cm respectively and for the 5pn Q n Rydberg state the 

results are 484.8cm 1 and 594.6cm 1 respectively. The difference of only about 111 

between theory and experiment for the lower two levels is quite encouraging 

considering that the matrix elements are quite small and only small valence space 

Cl wave functions have been used. Gal lusser and Dressier point out that the Q 

state matrix element is not as reliable as that for the C state which has the 

largest matrix element. The experimental Q state matrix element deviates by about 

3/2 . 2 

15Z from the value predicted from an n* behavior. For the L n repulsive state 
the calculated and experimental matrix elements with the C state are 546.7cm ^ and 
549.0cm * respectively while for the K state the results are 320cm * and 
-250+50cm ^ respectively and for the Q state we have 218.2cm ^ and (200cm ^). 
Gallusser and Dressier"* have indicated that the Q and K state results are not as 
reliable as the C state matrix elements with which we obtain excellent agreement. 


The NO calculations have been performed with a [3s,2p,ld] basis set for the 

valence orbitals supplemented by a set of 18 diffuse 2pJr^ Gaussian primitives 

(centered at the midpoint) for the Rydberg orbitals. The calculations give the 

contribution of only the 1=1 partial wave to the width matrix element. An 

2 

additional contribution to the Rydberg II ptate widths can be expected from the 1=2 

partial wave. The 1=2 contribution has been calculated by adding to the above 

basis a set of 12 3d Gaussians centered at the midpoint. The 1=2 functions 

xz 

increase the B state width by 3 .4% and the L state width by 1.7% • Therefore the 

2 

1=2 wave makes only a minor contribution to the NO Dwidths. 

IV Application to 


Using the above procedure for calculation of the width matrix elements 

3 

combined with previously calculated potential curves and the quantum defect theory 
of Giusti^. cross sections and rate constants have been determined for the DR of 
© 2 + leading to 0(^S), the upper state of the important atmospheric green line 
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emission. The calculations include the effect of autoionization, i.e. ejection of 
the electron after capture and before dissociation occurs. The calculations do not 
include the effect of indirect recombination through intermediate Rydberg states. 
Indirect recombination is expected to be important for 0(*S) generation and is the 
subject of current work. 

Fig. 1 shows the calculated rates for particular ion vibrational levels. 
Clearly, the rate is very sensitive to ion vibrational population. The rates for 
v=2, 3, and 5 are more than two orders of magnitude greater than the rate for v=0 
at low electron temperatures . 

Fig. 2 shows the calculated rates as a function of vibrational temperature. 
As expected from Fig. 1, the rates are sensitive to vibrational temperature at low 
electron energies but insensitive at high electron energies where Franck-Condon 
factors become small and autoionizaticn is important. Note that the electron 
temperature at which the rates become insensitive to vibrational temperature is 
highly dependent on the shapes of the potential energy curves. 

DR calculations over a wide range of temperatures for the remaining 
dissociative routes of 0 ^ are currently in progress. 
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